Hidden glassy behaviors in an ideal Heisenberg Kagome antiferromagnet 
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Dynamics of classical Heisenberg spins, Si = (si,Si Z ), on the Kagome lattice has been studied. 
An ideal Heisenberg Kagome antiferromagnet is known to remain disordered down to T — due to 
the macroscopic degeneracy of the ground state. Through the study, however, we find that Si and 
their planar components Sj behave in a qualitatively different way and especially planar spins (sj) 
show the exotic glass-like transition in the very low temperature T ~ 0.003 J (J: spin exchange). 
The glassy behaviors of s* would be found to be driven by the spin-nematic fluctuations, different 
from ordinary spin glasses by disorders or anisotropies. 



Magnetic systems in which the exchange energy is not 
minimized by a simple regular arrangement of the spins 
are called the frustrated magnetic systems. Frustrated 
antiferromegnets induced by the geometry of the lattice 
have been a recent subject attracting much experimental 
and theoretical interests 0. One example of such sys- 
tems is the antiferromagnet on the Kagome lattice. It 
consists of triangles. The Kagome lattice is, however, 
more frustrated than the triangular lattice because tri- 
angles on the former lattice share only one vertex, but 
on the latter share one side. Indeed, classical Heisen- 
berg spins on the Kagome lattice have an infinitely de- 
generate ground state and some novel ground state is 
expected as T — > without developing long-range order 
at any T, which has been a main reason drawing consid- 
erable attentionj^OliELIHEi- Their experimental realiza- 
tions are, for instance, the insulating layered compound 
SrCr 8 _ x Ga4+xOi 9 (Cr 3+ carries S = 3/2) 0, HQ and 
the jarosites AFe 3 (OH) 6 (S0 4 )2 (A=K_H 3 0, Na, NH 4 , 
and so forth; Fe 3 + carries S = 5/2) Eft [HI- 

Extremely high degeneracy of the ground state inhibits 
the long-range order at any T. But it has been known 
that thermal fluctuations in the system select a subset of 
ground state manifold at low T (known as ordering by 
disorder) 0, |(| , even if this subset is not clear enough to 
lead to long-range Neel order. Chalker et a/.jj| have ar- 
gued by low temperature expansions that thermal fluctu- 
ations resolve degeneracy of ground states and cause the 
system to select a coplanar nematic ground states. There 
are at least two coplanar ground states possessing long- 
range Neel order; one is the so-called q — state and the 
other is \/~3 x y3 state (See the inset in Fig^a); in the 
q = state, spins along lines are placed in an alternating 
sequence compared to the v3 x v3 state) . Reimers and 
Berlinsky[^ have shown, through the Monte Carlo sim- 
ulation, that an increase in the order parameter of the 
V3 x V3 state is dominant over the q = state. In addi- 
tion, Sachdev^3 has also argued that the \fi x \fi state 
is further accepted and favored in the low temperature 
limit by quantum fluctuations. 

Magnetic frustration used to be found with the site 
disorder, even if these two features may appear indepen- 
dently. When both, disorder and frustration, are strong, 
its interplay gives place to the spin glass phase^S]- On 
the other hand, the antiferromagnetic Heisenberg model 



on the Kagome lattice is an ideal geometrically frus- 
trated magnet without disorders. Nevertheless, in some 
of prototype systems of the Kagome lattice, the spin 
glass ordering at a finite temperature is actually ob- 
served because of deviations from an ideal one. In 
SrCr 8 _ :! .Ga4 +a ;Oi9, T g = 3.3 K is found in a case of x = 
corresponding to 89% occupation of Cr sites 0, @, 0- 
In the system, Cr 3+ just enters three different octa- 
hedral sites (namely, 2a, 12k, 4f) with an almost ran- 
dom distribution and SrCr 8 Ga 4 0i9 has then ~14% site 
disorder (nonmagnetic impurities) on the Kagome layer 
(12k sites). The glass order is enhanced by the ad- 
ditional dilution of Cr sites 0- The glassy phase in 
nondisordered frustrated magnetic systems has been at- 
tracting much interests. Most two-dimensional mag- 
netic materials exhibit some kinds of spin anisotropy. A 
critical slowing down of magnetic fluctuations is found 
in the hydronium-jarosite (H30)Fe3(OH)g(S04)2 with 
much less disorders 11]. The glass phase in the sys- 
tem may be a bit different from in SrCrgGa40ig, due to 
the xy anisotropy ("easy-plane"), which could drive the 
Kagome antiferromagnet to experience the topological 
Kosterlitz-Thouless-type transition to a glassy stated- 
The glassy behavior in the Kagome antiferromagnet with 
the easy-axis type anisotropy has been also explored p^]. 

In this paper, through the study of dynamical proper- 
ties, we report the hidden glassy behaviors of Heisenberg 
spins on the ideal Ka gom e lattice without any disorder 
or anisotropy. Keren|15j has investigated the stability 
of the spin order against excitation and the dynamical 
responses of spins on the Kagome lattice. Actually, in 
his study, there is found no transition down to T = as 
there should not be. Nevertheless, we find that behaviors 
of Si and its planar component s, are qualitatively differ- 
ent from each other and Sj undergoes a glassy transition 
around T ~ 0.003 J unlike S,-. 

Taking the classical Heisenberg spins as Si = (s,, Si Z ), 
we have the Hamiltonian 



H = J^Si ■ s 3 + J^2s iz S 



(1) 



where sf + Sf z = 1. Reimers and Berlinsky 6] have 
demonstrated that the tendency toward spin-nematic or- 
der starts at T ~ 0.01J, which motivates us to explore 
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FIG. 1: Local spin-spin correlation functions at ui = with re- 
spect to temperatures T for both cases of relaxing and anneal- 
ing, (a) gives behaviors of J(0) which is history-independent, 
whereas (b) shows that j(0) is history-dependent approxi- 
mately below T ~ 0.003J. The inset in (a) is the spin con- 
figuration for the \/3 x \/3 state and the inset in (b) is the 
profile of J(0) — j(0), i.e. J z (0) with T (the line is a guide for 
the eye) . 



the respective dynamics of Si and at the relevant tem- 
perature range, i.e. T < 0.01J. Particularly, we have 
interests in the local spin-spin correlation functions as 
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The correlation functions such as J {to) and j(u>) 
can be most properly evaluated by the Monte Carlo 
simulation 1 16|. The time evolution of spins are governed 
by the equation of motion dSi/dr — — JS.; x Yjj& Sj< A 
practical integration of the equation is done using the 
Suzuki- Trotter decomposition on the three sublattices 
(say A, B, and C) of the Kagome lattice, {Si(r + 5t)} = 
e^ A+B+c ^ T {S l (T)}, where e ^+B+c)Sr is decomposed 



(X+Y)St 



n 



i=l 



using e 

and pi = p? = j»4 = P5 = = 1/(4 — 4 1 / 3 ) and 
pa = 1 — 4p[lq. We use the hybrid Monte Carlo proce- 
dure which combines the Metropolis update and the over- 
relaxation update^! • One hybrid Monte Carlo step in 
our calculation consists of two Metropolis steps and four 
overrelaxation steps. Using the hybrid algorithm, O(10 4 ) 
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FIG. 2: Comparison of inverse intensities 1/J(0) and l/j(0) 
with temperatures T. Inverse intensities are well fitted with 
a linear T. l/j(0) shows a clear refraction near ~ 0.003 J, 
which is not in 1/J(0). 



hybrid steps are typically used to reach the equilibrium 
configuration^^. A system size is defined by the linear 
dimension L of the lattice and the total number of spins 
N is 3L 2 . In an actual simulation for dynamical correla- 
tion functions, the time integration has been performed 
up to T max = 26214. 4J _1 with a step of St = 0.2J -1 . 
The thermal average of time-dependent variables is com- 
puted over the observables obtained by evolving all inde- 
pendent initial equilibrium configurations. Typically, 10 3 
configurations were generated at a given T. An adopted 
system size for dynamical calculations is L = 24, that is, 
1728 spins[H. 

In Fig^ the temperature-dependences of J(0) and j(0) 
are provided both when relaxing and annealing the sys- 
tem. J(0) (or j(Q)) is a quasi-elastic response of local 
spin-spin correlation function and a quantity directly re- 
lated to the spin-lattice relaxation rate T\ of a local probe 
like the nuclear magnetic resonance (NMR). For a case 
of relaxing, the nominal state at T = is taken as the 
V3 x v3 state. For a case of annealing, on the other 
hand, after heating the y/3 x -\/3 state to T = 0.05 J, we 
obtain the equilibrium state by gradual cooling. Fig^] 
shows that J(0) is history- independent, whereas j(0) is 
history-dependent roughly below T ~ 0.003J. History- 
dependent magnetic responses can be one of indications 
of the glassy phase because of the slow dynamics of 
spins in the phase. Hereafter the aimed configuration 
is obtained by warming the system from the x \fZ 
state, if not stated otherwise. In the inset of the figure, 
J(0) — i(0), i.e. J z (0) are given. A much sharper com- 
parison between J(0) and j(0) is done by taking their 
inverses. In FigO 1/J(0) and l/j(0) are shown more or 
less linear with respect to T and l/j(0) suffers an abrupt 
and appreciable change of the proportionality constant 
near T ~ 0.003 J, which is not found in 1/J(0). 

It should be also interesting to investigate the behav- 
iors of J(lo) and j(u>) in the low energy region (lj <C 
0.003 J). Fischer and Kinzel 20] have extended the dy- 
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FIG. 3: J(w)/J(0) and j(w)/j(0) for small u) 's (|w| < 0.003J) 
at several temperatures. In j(u))/j(0), qualitative changes of 
w-behaviors are observed around T ~ 0.003J. 



namical analysis of early spin glass theories and pre- 
dicted a crossover in low energy behaviors of the ac mag- 
netic susceptibility (x"(w); basically, same as J(u>)) from 



AT) 



near 



x"(w) ~ uj in high temperatures to x"( w ) 
T 9 and especially ^(T 9 ) = 1/2. Paulsen aL[21| have 
found that x"( w ) m the spin glass system Euo.4Sro.6S has 
a power-law dependence on frequency like x"( w ) ~ w^ T ' 
and f(T) has abruptly changed from v(T < T s ) ~ 0.09 
to v(T > T s ) ~ 0.4 around the freezing temperature T g . 
Broholm et a£.|S| have also reported the similar power 
law behavior like u"' T ' of the local response function in 
SrCr 8 _ 2 ,Ga4 +:E Oig and changes of v(T) around T g , that 
is, y{T) ~ below T g and ~ 1 in high temperatures. In 
FigEl the low energy responses of the spin-spin correla- 
tion functions are given. Behaviors of j(u>) in the lower 



panel of Fig is appealing in that 



expected to follow - 1/(uj 2 + T(T) 2 )) below T ~ 0.003 J 
changes to ~ u v (v > 0) above T ~ 0.003 J. No 
similar changes are found in J(w) as shown in the upper 
panel of the figure. All previous results of Figs lll3l on 
j(oj) are subordinate to the glass-like phase transition at 
T ~ 0.003 J and quite consistent with each other. On the 
contrary, any indication of spin phase transition down to 
T — > is not found in magnetic responses of SL (i.e. J(0) 
and J(u>)), which is consistent with Keren'sfl^. 

We now importantly note that Si and Sj behave in a 
qualitatively different way. In order to understand the 
differences, we pay attention to the low energy sector 
of the Hamiltonian. Following Chalker et aL(4|, Ti, in 
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FIG. 4: Figurative diagram of the Hamiltonian (TL) in the 
low-temperature expansion and its planar counterpart (TL s i 
and TL s ii). TL s i and TL s n , projected to the plane of = 0, are 
for two different local spin distortions, while TL has a stable 
saddle point for (small) spin distortions. 



the low temperature expansion, on the Kagome lattice is 
schematically 

H = Hq+ Wh(4l e p) + "ftanli(ea£/3, e /3, ■■■), 

where e a (ep) is the fluctuation perpendicular to Si in 
the particular planar ground state and in (out of) the 
spin plane and Hh and H an h includes harmonic and an- 
harmonic fluctuations, respectively. In harmonic fluctua- 
tions (within the linear spin waves), e 2 a gives three trans- 
verse modes and e| two transverse modes and one zero 
mode per unit cell comprising three spins. The role of 
quartic potential for the zero mode is then unique and 
is to stabilize the planar states by inducing the spin- 
nematic fluctuation. Now, let us think of the planar 
counterpart of TL, i.e. TL S , 

In TL S , terms of e Q e| and disturb the system like disor- 
ders (but not static) in a complicated way. For a weakly 
(or non-) disordered case, it is known that the ground 
state configuration is such that the energy of each sepa- 
rate triangle is minimized ("rule of satisfied triangles"), 
implying the glass phase will not occur |22|. In Ti s , how- 
ever, in Fig0J the energy minimization of each triangle 
depends on the local spin distortions produced by neigh- 
boring triangles, in other words, by correlations between 

■2 ■ 2 ■ 2 

fluctuating degrees of freedom (e.g. e^Cg or ). Al- 

ternatively, terms of e a e 2 p and ejg can be also understood 
to lead to planar spins of dynamically varied magnitudes, 
which act as strong disorders such that glassy behavior 
may occur. To integrate out {e a } in the partition func- 



tion or correlation functions leaves terms of e^. If we 
note its role in TL and the ground state, the glass-like or- 
der of planar spins can be attributed to the spin-nematic 
fluctuations. The glass-like order in the system is differ- 
ent from conventional spin glasses in that it is originated 
by the dynamical fluctuations of inherent degrees of free- 
dom, not by the static physical disorders or anisotropics 
in the systems. 
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FIG. 5: Static spin susceptibilities for Si and Si, xs and Xs- 
The inset gives the spin susceptibility for XY spins, xxy- 
Calculations are done for 2700 spins on the Kagome lattice 
(L = 30). 

Finally, it may be interesting to test a case of 
XY spins (planar spins by definition) without such 
fluctuations 23] , for which the easiest way is probably 
to examine the static spin susceptibility. It is also a 
good indicator of the spin ordering and shows a cusp 
at the transition temperature in the spin glass 124) . 
The spin susceptibilities are defined for Si and Sj as 
XS (T) « [{Si ■ Sj) - (Si) ■ (Sj)] and Xs (T) ex 

1 / T J2i^j i( s i ■ s j> - ( s i) ■ ( s j)]. respectively. In FigEJ it 
is found that only Xs(T) shows a cusp at T ~ 0.003J 



and other xs (T) and xxy (T) increases smoothly without 
any singularities with respect to T, even if all the three 
cases are started from the seemingly same spin configu- 
ration (i.e. V3 x V3 state). Spin susceptibilities in Fig|3 
ascertain the glassy behavior of the planar components 
of Heisenberg spins on the Kagome lattice and further 
support that it would be induced by the spin-nematic 
fluctuation. As mentioned, one of the most outstanding 
features of the frustrated magnet is the high degeneracy 
of the ground state, whose peculiarity could result in the 
spin-nematic ordering in the limit of very low tempera- 
tures (T <C J). The nematic glassy state should then 
be crucial and eminently important in understanding the 
novel ground state of the frustrated magnet. 

In summary, we have found and discussed qualitative 
differences between S; and its planar component Sj in an 
ideal Heisenberg Kagome antiferromagnet. It is remark- 
able to find the hidden glass-like order of Sj even in the 
Heisenberg model on the Kagome lattice without any dis- 
order or anisotropy. The glassy behavior has been shown 
through dynamical spin responses and static susceptibil- 
ities. We find that it would be driven by the dynamical 
fluctuation of spin-nematic order, different from ordinary 
spin glasses by static disorders or anisotropies. 
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ported by the U.S. Department of Energy under contract 
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